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We report results from a search for production of a neutral Higgs boson in association with a b 
quark. We search for Higgs decays to r pairs with one r subsequently decaying to a muon and the 
other to hadrons. The data correspond to 2.7 fb _1 of pp collisions recorded by the DO detector at 
^J~s~ — 1.96 TeV. The data are found to be consistent with background predictions. The result allows 
us to exclude a significant region of parameter space of the minimal supersymmetric model. 

PACS numbers: 12.10.Dm, 12.60.Pr, 12.60.Jv 



The current model of physics at high energies, the 
standard model (SM), has withstood increasingly pre- 
cise experimental tests, although the Higgs boson needed 
to mediate the breaking of electroweak symmetry has 
not been found. Despite the success of the SM, it has 
several shortcomings. Theories invoking a new fermion- 
boson symmetry, called supersymmetry [l| (SUSY), pro- 



vide an attractive means to address some of these includ- 
ing the hierarchy problem and nonunification of couplings 
at high energy. In addition to new SUSY-specific part- 
ners to SM particles, these theories have an extended 
Higgs sector. In the minimal supersymmetric standard 
model (MSSM) there are two Higgs doublet fields which 
result in five physical Higgs bosons: two neutral scalars 
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(h, H), a neutral pseudoscalar (A) and two charged Higgs 
bosons (H ). The mass spectrum of the Higgs bosons 
is determined at tree level by two parameters, typically 
chosen to be tan (3, the ratio of the vacuum expectation 
values of up- type and down- type scalar fields and Ma, 
the mass of the physical pseudoscalar. Higher order cor- 
rections are dominated by the Higgsino mass parameter 
fi and the mixing of scalar top quarks. 

In this Letter, we present a search for neutral Higgs 
bosons (collectively denoted </>) produced in association 
with a b quark. The specific Higgs boson decay mode 
used in this search is </>—»■ rr with one of the r leptons 
subsequently decaying via t — > fifi-v^ (denoted t^) and 
the second via r — > hadrons + v r (denoted t^). In the 
MSSM the Higgs coupling to down-type fermions is en- 
hanced by a factor oc tan (3 and thus the Higgs production 
cross section is enhanced by a factor oc tan 2 f3 relative to 
the SM, giving potentially detectable rates at the Teva- 
tron. Two of the three neutral Higgs bosons have nearly 
degenerate masses over much of the parameter space, ef- 
fectively giving another factor of two in production rate. 
A previous search in this final state was carried out by 
the DO experiment Searches in the complementary 



rithm [&] with a radius of R 
R 



0.5 in (y, ip) space, where 



bbrr, rrbb Q , 



[4, 5], and 



channels 

<f)b — + bbb [6, 7] have also been carried out by the LEP, DO, 
and CDF experiments. By searching in complementary 
channels we reduce overall sensitivity to the particular 
details of the model. The brr final state is less sensitive 
to SUSY radiative corrections than the bbb final state, 
and has greater sensitivity at low Higgs mass than the 
<f) — > tt channel, as the &-jet in the final state reduces the 
Z* — > tt background. Furthermore, an additional com- 
plementary channel will contribute to an even stronger 
exclusion when combining different searches. The result 
presented in this Letter uses an integrated luminosity of 
2.7 fb _1 which is eight times that used for the previous 
result in this channel. Because of analysis improvements, 
the gain in sensitivity compared to the prior result is 
greater than expected from the increased integrated lu- 
minosity only. We also extend the Higgs mass search 
range relative to the previous result in this channel. 

The DO detector [8] is a general purpose detector lo- 
cated at Fermilab's Tevatron pp collider. The Tevatron 
operates at a center of mass energy of 1.96 TeV. This 
analysis relies on all aspects of the detector: tracking, 
calorimetry, muon detection, the ability to identify de- 
tached vertices and the luminosity measurement. 

This search requires reconstruction of muons, hadronic 
r decays, jets (arising from b quarks) and neutrinos. 
Muons are identified using track segments in the muon 
system and are required to have a track reconstructed 
in the inner tracking system which is close to the muon- 
system track segment in r\ and ip. Here r\ is the pseu- 
dorapidity and <p is the azimuthal angle in the plane 
perpendicular to the beam. Jets are reconstructed from 
calorimeter information using the DO Run II cone algo- 



y/ (Ay 2 + Ap 2 ) and y is the rapidity. Jets are ad- 
ditionally identified as being consistent with decay of a 
fo-flavored hadron (6-tagged) if the tracks aligned with 
the calorimeter jet have high impact parameter or form 
a vertex separated from the primary interaction point in 
the plane transverse to the beam as determined by a neu- 
ral network (NNf,) algorithm [l(|. Hadronic r decays are 
identified [1 ill as clusters of energy in the calorimeter re- 
constructed [9J using a cone algorithm of radius R = 0.3 
which have associated tracks. The r candidates are then 
categorized as being one of three types which correspond 
roughly to one-prong r decay with no 7r°s (called Type 1), 
one-prong decay with 7r°s (Type 2) and multiprong de- 
cay (Type 3). A final identification requirement is based 
on the output value of a neural network (NN T ) designed 
to separate r leptons from jets. The missing transverse 
energy 1£t is used to infer the presence of neutrinos. The 
$t is the negative of the vector sum of the transverse 
energy of calorimeter cells satisfying \r)\ < 3.2. J^t is cor- 
rected for the energy scales of reconstructed final state 
objects, including muons. 

Signal acceptance and efficiency are modeled using 
simulated SM <j)b events generated with the PYTHIA event 
generator [l2T ] requiring the b quark to satis fy p t > 
15 GeV/c and 1 77 1 < 2.5 and using the CTEQ6L1 JlJ par- 
ton distribution functions (PDF). The TAUO LA Jl4l| pro- 
gram is used to model r decay and EVTGEN [l5| is used 
to decay b hadrons. The dependence of the Higgs boson 
decay width on tan (3 is included by reweighting pythia 
samples, and the kinematic properties are reweighted to 
the prediction of the NLO program mcfm 16]. The 
generator outputs are passed through a detailed detec- 
tor simulation based on GEANT [lTj. Each GEANT event 
is combined with collider data events recorded during a 
random beam crossing to model the effects of detector 
noise, pileup, and additional pp interactions. The com- 
bined output is then passed to the DO event reconstruc- 
tion program. Simulated signal samples are generated for 
different Higgs masses ranging from 90 to 320 GeV/c 2 . 

Backgrounds to this search are dominated by Z+jets, 
tt, and multijet (MJ) production. In the MJ background 
the apparent leptons primarily come from semilcptonic 
b hadron decays, not r decays. Additional backgrounds 
include H^+jets events, SM diboson production and sin- 
gle top quark production. Except for the MJ contribu- 
tion, all background yields are estimated using simulated 
events, with the same processing chain used for signal 
events. The Z+.jets, IF+jets and tt samples are gener- 
ated using ALPGEN [10| with pythia used for fragmenta- 
tion. The diboson samples are generated using PYTHIA. 
For simulated samples in which there is only one lepton 
arising from the decay of a W boson or from tt — » ^+jets, 
the second lepton is either a jet misidentified as a r or a 
muon+jet system from heavy flavor decay in which the 
muon is misidentified as being isolated from other activ- 
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ity. 

Corrections accounting for differences between data 
and the simulation are applied to the simulated events. 
The corrections are derived from control data samples 
and applied to object identification efficiencies, trig- 
ger efficiencies, primary pp interaction position (primary 
vertex) and the transverse momentum spectrum of Z 
bosons. After applying all corrections, the yields for sig- 
nal and each background are calculated as the product 
of the acceptance times efficiency determined from sim- 
ulation, luminosity and predicted cross sections. 

The initial analysis step is selection of events recorded 
by at least one trigger from a set of single muon triggers 
for data taken before the summer of 2006. For data taken 
after summer 2006 we require at least one trigger from 
a set of single muon triggers and muon plus hadronic r 
triggers. The average trigger efficiency for signal events 
is approximately 65% for both data epochs. 

After making the trigger requirements a background- 
dominated pre-tag sample is selected by requiring a re- 
constructed primary vertex for the event with at least 
three tracks, exactly one reconstructed hadronic r, ex- 
actly one isolated muon, and at least one jet. This anal- 
ysis requires the r candidates to satisfy Ej, > 10 GeV, 
p T T > 7(5) GeV/c and NN T > 0.9 for Type 1(2) taus, 
E T T > 15 GeV, p T T > 10 GeV/c and NN T > 0.95 for 
Type 3 taus. Here E7p is the transverse energy of the 
r measured in the calorimeter, pTp is the transverse mo- 
mentum sum of the associated track(s). The muon must 
satisfy pl^ > 12 GeV/c and \t)\ < 2.0. It is also required 
to be isolated from activity in the tracker and calorime- 
ter 0. Selected jets have E T > 15 GeV, \r)\ < 2.5. The 
t, the muon and jets must all be consistent with aris- 
ing from the same primary vertex and be separated from 
each other by AR > 0.5. In addition, the muon and t are 
required to have opposite charge, and the (p,$t) mass 

variable M = y^r^/p^ 1 ~ cos ( a ¥>(m,$t))) must 
satisfy M < 80, 80, and 60 GeV/c 2 for events with rs 
of Type 1, 2 and 3 respectively. Here E^ is the energy of 
the muon, and A.tp is the opening angle between the $t 
and muon in the plane transverse to the beam direction. 

A more restrictive &-tag subsample with improved sig- 
nal to background ratio is defined by demanding that at 
least one jet in each event is consistent with b quark pro- 
duction [10|| . The 6-jet identification efficiency in signal 
events is about 35% and the probability to misidentify a 
light jet as a b jet is 0.5%. 

All backgrounds except MJ are derived from simulated 
events as described earlier. The MJ background is de- 
rived from control data samples. A parent MJ-enriched 
control sample is created by requiring a muon, r, and 
jet as above, but with the muon isolation requirement 
removed and with a lower quality (0.3 < NN T < 0.9) r 
selected. This is then used to create a 6-tag subsample 
which requires at least one of the jets to be identified as a 



b jet with the same b jet selection as earlier. The residual 
contributions from SM backgrounds are subtracted from 
the MJ control samples using simulated events. 

To determine the MJ contribution in the pre-tag anal- 
ysis sample, a data sample is used that has the same 
selection as the pre-tag analysis sample except that the 
muon and r charges have the same sign. This same-sign 
(SS) sample is dominated by MJ events. After making 
a subtraction of other SM background processes which 
contribute to this sample, the number of MJ events in 
the opposite-sign (OS) signal region is computed by mul- 
tiplying the SS sample by the OS/SS ratio, 1.05 ± 0.02, 
determined in a control sample selected by requiring a 
non-isolated muon.. 

For the 6-tag analysis sample, statistical limitations 
require a different approach for the MJ background eval- 
uation than for the pre-tag sample. For the 6-tag sample, 
two methods are used. For the first method, the per jet 
probability P tag that a jet in the SS MJ control sub- 
sample would be identified as a & jet is determined as a 
function of jet px- We apply P tag to the jets in the SS 
pre-tag sample to determine the yield in the &-tag sam- 
ple. For the second method, the MJ background is deter- 
mined by multiplying the 6-tag MJ control sample yield 
by two factors: (1) the probability that the non-isolated 
muon would be identified as isolated, and (2) the ratio of 
events with a r candidate passing the NN T requirements 
to events with r candidates having 0.3 < NN T < 0.9 
as determined in a separate control sample. The final 
MJ contribution in the 6-tag analysis sample is deter- 
mined using the MJ shape from the first method with 
the normalization equal to the average of the two meth- 
ods. We include the normalization difference between 
the two methods in the systematic uncertainty on the 
MJ contribution. 

The signal to background ratio is further improved us- 
ing multivariate techniques. Two separate methods are 
used, one to address the tt background and one to re- 
duce the MJ background. For the tt background, a neu- 
ral network (NNt op ) is constructed using Ht = Y>j e tsET, 
E to t = ^jetsE+ET+E^, the number of jets and A<p(/i, r) 
as inputs. For the MJ background, a simple joint likeli- 
hood discriminant [LLmj) is constructed using p T , p T , 
AJ?(/x,r), M MT and M„ T!/ . Here M„ T denotes the invari- 
ant mass of the muon and tau, and M„ TI/ is the invariant 
mass computed from the muon, t, and $t momentum 
vectors. The final analysis sample is defined by select- 
ing rectangular regions in the NN top versus LLmj plane. 
The regions have been identified for each r type and 
each Higgs boson mass point separately by optimizing 
the search sensitivity using simulated events. The signal 
to background ratio improves by up to a factor of two 
when applying these requirements. 

Table U shows the predicted background and observed 
data yields in the analysis samples. Between 5% and 10% 
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Pre-tag 


6-tagged 


Final 


tt 


66.0 ± 1.3 


39.6 ± 0.8 


20.3 ± 0.6 


Multijet 


549 ± 26 


38.5 ± 2.3 


28.1 ± 1.9 


Z(— » tt) + jets 


1241 ±8 


18.8 ±0.3 


16.3 ±0.3 


Other Bkg 


267 ±6 


5.1 ±0.1 


4.1 ± 0.1 


Total Bkg 


2123 ±28 


102 ±2.4 


68.8 ±2.0 


Data 


2077 


112 


79 


Signal 


14.4 ±0.3 


4.8 ±0.1 


4.6 ±0.1 



TABLE I: Predicted background yield, observed data yield 
and predicted signal yield and their statistical uncertainties 
at three stages of the analysis. The signal yields are calculated 
assuming tan/3 = 40 and a Higgs mass of 120 GeV/c 2 for the 
m^ ax and n = -200 GeV/c 2 



scenario. 



of <j) — » Tfj,Th decays are selected depending on M^. 

Systematic uncertainties arise from a variety of 
sources. Most are evaluated using comparisons between 
data control samples and predictions from simulation. 
The uncertainties are divided into two categories: (1) 
those which affect only normalization, and (2) those 
which also affect the shape of distributions. The sources 
in the first category include the luminosity (6.1%), muon 
identification efficiency (4.5%), tk identification (5%, 4%, 
8%), Th energy calibration (3%), the it and single top 
cross sections (11% and 12%), diboson cross sections 
(6%), Z+(u,d,s,c) rate (±2%, -5%) and the W + b and 
Z + b cross sections (30%); those in the second include jet 
energy calibration (2%-4%), 6-tagging (3%-5%), trigger 
(3%-5%), and MJ background (33%, 12%, 11%). For 
sources with three values, the values correspond to r 
Types 1, 2 and 3 respectively. 

After making the final selection, the discriminant D 
is formed from the product of the NN top and LLmj 
variables, D = LLmj x NN top . The resulting distribu- 
tions for the predicted background, signal and data are 
shown in Fig. [lja). This distribution is used as input 
to a significance calculation using a modified frequentist 
approach with a Poisson log-likelihood ratio test statis- 
tic [20j • In the absence of a significant signal we set 95% 
confidence level limits on the presence of neutral Higgs 
bosons in our data sample. The cross section limits are 
shown in Fig. [ljb) as a function of Higgs boson mass. 
These are translated into the tan (3 versus Ma plane in 
the m% ax , n = -200 GeV/c 2 MSSM benchmark scenario 
(2lT |. giving the excluded region shown in Fig. QJc). The 
signal cross sections and branching fractions arc com- 
puted using feynhiggs [II]. Instabilities in the theo- 
retical calculation for tan (3 > 100 limit the usable mass 
range in the translation into the (tan/3, Ma) plane. 

In summary, this Letter reports a search for produc- 
tion of Higgs bosons in association with a b quark using 
eight times more data than previous results for this chan- 
nel. The data are consistent with predictions from known 
physics sources and limits are set on the neutral Higgs 



boson associated production cross section. These cross 
section limits, a factor of three improvement over previ- 
ous results, are also translated into limits in the SUSY 
parameter space. 
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